1 Introduction 



With the advent of very energetic hadronic colliders like LHC and SSC, 
whose aim is to have access to the mechanism of electroweak symmetry break- 
ing, a precise estimate of inclusive production of neutral clusters is mandatory 
to pin down signals due to Higgs particle or new physics Perturbative QCD 
is the appropriate framework to perform such a calculation provided large trans- 
■ fer momenta are involved. For this purpose leading order (hereafter denoted 

Q\ • as LO ) predictions - based on evaluations of partonic cross sections at tree 

level and evolution of structure and fragmentation functions at one loop level 
- are too rough. A consistent calculation at next to leading order (hereafter 
denoted as NLO) needs a NLO evaluation of parton-parton subprocesses which 
has been performed by our group a few years ago utilizing existing results 
on the O(ajj) matrix element [|| and two loop evolved structure and fragmen- 
tation functions. Various sets of structure functions based on a NLO analysis 
using deep inelastic lepton nucleon scattering data, Drell-Yan production and 
J> ' direct photon production at hadronic accelerators are available in the literature 

ill- 

(N 

Up to now such an analysis has not been performed for pions fragmentation 
functions since the only available derivation at NLO has been done for heavy 
quarks |7|. The parametrizations presently available || are based on a LO 
analysis of rather old e + e~ data and semi inclusive deep inelastic muon nucleon 
Q_i| scattering. Our aim is to perform a complete NLO evaluation of neutral pions 



inclusive production from hadronic colliders in order to estimate, as precisely as 



possible, the tt° rates at LHC and SSC. The first step will consist in performing 
an extraction of pions fragmentation functions at NLO using e + e~ data and 
hadronic data on one particle inclusive production. This will be done using 
three different methods. The first one, which does not exactly correspond to 
a NLO analysis but rather to an improved LO approximation, is based on the 
Monte-Carlo simulator HERWIG Q, whereas the second and the third ones are 
obtained from a two loop evaluation of evolution kernels previously computed 



1C , Ojtogether with NLO calculation of one hadron inclusive production from 
e + e~ |12| and hadronic colliders M using respectively natural and optimized 
scales. As we will see it is not possible to extract an unique set of fragmentation 
functions which fits e + e~ data around y/S = 30 GeV |l5[ 16, 17, 18] and 



hadronic data from fixed target domain |T^, to collider range ® n m. 
We will therefore take different sets corresponding to different hypotheses on 
input fragmentation functions. 



The paper is organized as follows. We recall the expressions of one parti- 
cle inclusive production from e + e~ and hadronic collisions and also give the 
evolution equations for fragmentation functions in section 2. Then in section 3 
we discuss the extraction of various sets of fragmentation functions, first from 
HERWIG simulation and after through an exact NLO derivation. Predictions 
at LHC using the different sets previously obtained are displayed in section 4 
together with a discussion of the resulting theoretical uncertainty. We give our 
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conclusions in section 5. 



One particle inclusive production at next-to-leading 
order 



Let us consider the inclusive production of a hadron H via the generic 
reaction A + B — » H where A and B stand for hadrons and/or leptons. The 
cross-section can be written as a convolution of the fragmentation functions 
Df 1 (z, M 2 ) with the partonic cross-section: 
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where zh is the reduced energy of the hadron H: zh = lEuj^fS and is the 
scattered angle of the parton 1. The inclusive production of the parton 1 via the 
reaction A + B —* I has the following perturbative development: 
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Finally Df 1 (z, M 2 ) represents the number of hadrons H inside the parton 1 
carrying the fraction of impulsion z from H, evolved at the scale M 2 . These 
fragmentation functions satisfy Altarelli-Parisi type evolution equations: 
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The evolution kernels have the perturbative development: 
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In the following, we will drop the superscript index T. The Altarelli-Parisi 
kernels have been computed up to two loops order by Curci, Furmanski and 
Petronzio [11]. In the LO approximation one keeps only the first order in 
the perturbative development of the partonic cross-section and in the evolution 
kernels whereas at NLO one keeps the first and second terms in the perturbative 
expansion for both partonic cross-section and evolution kernels. We can split 
these fragmentation functions into a non-singlet and a singlet part: 
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Dt(z,M}) = 2 
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D s (z,M 2 f ) = ^2(D*(z,M 2 f )+D?(z,M 2 f 
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In the evolution equations the singlet part D$ is coupled to the gluon fragmen- 
tation function whereas the non-singlet parts D~ and D + are decoupled. Note 
a misprint in ref. [jlTj] and the correspondence with our notation: 
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Once input fragmentation functions have been specified at some reference scale 
Mjo the evolution equations are solved using an inverse Mellin transform tech- 
nique. Let us consider now in detail the partonic cross-sections. 



2.1 e + e 



7T 



The partonic cross sections from e + e collisions read at next-to-leading 
order: 
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where <ro is the usual point like cross-section 

47ra 2 



(9) 



0Q 



3Q 2 



a is the QED coupling constant and Q 2 is the invariant mass of the e + e pair. 
The functions KT , 
|| (see also @). 



The functions KT, , iCj and Kg have been extracted from the reference 
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In the above equations two scales are involved: the renormalization scale \x 
at which the running coupling constant a s is evaluated and the fragmentation 
scale Mf at which fragmentation functions are evolved. The choice for these 
scales is rather arbitrary. Note that for every y, Kj(y) is negative, so the 
choice M 2 = Q 2 leads to a negative contribution to the partonic cross-section 

Eg da e+e -^ g /d 3 Pg. 

The running coupling of QCD a s is defined at the next-to-leading logarithm 
approximation by the approximate analytical formula: 
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In section 3 we will also use for a s the numerical solution of the equation: 
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which is more appropriate than eq.(14) for small scales fi. Indeed for large [i 
the two definitions agree but for small /x they can differ by more than 20 %. 
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The partonic cross-sections for hadronic collisions are given by 
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The variables V, W are defined by 
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and we also have 

VW l-V 



x\ = , x 2 



vw 1 — V 

and s = X1X2S. At NLO sixteen subprocesses contribute to the cross-section. 
The terms a correspond to the lowest order 2^2 parton scattering subpro- 
cesses whereas the terms K contain the one loop corrections to these subpro- 
cesses. In the hadronic case, we have three scales: the renormalization scale 
/i, the factorization scale for the initial state M ( the scale of the distribution 
functions) and the factorization scale for final state Mf (the scale of the frag- 
mentation functions). Schematically, the hadronic cross-section can be written 
as: 
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We show explicitly the dependence of the hadronic cross-section upon the three 
scales fx, M and M*. The four functions A, B, C and D depend on the scales M 
and Mf via the structure and fragmentation functions. In addition, A, B and C 
are scheme independent. We always use the MS scheme for final factorization 
whereas the initial factorization scheme is fixed by the set of structure functions 
used. 

Let us discuss now the partonic cross-sections. In order to determine the 
kinematical region where each partonic reaction dominates we have plotted 
in figures la, lb, lc and Id the partonic cross-sections E[ da p+p ->i/d 3 Pi for 
/ = g, u + u + d + d, s + s + c + c against Pt at the leading log level for 
various center of mass energies. We think it is meaningless to use next-to- 
leading formulae since the dependence on ln(M?) is not balanced. We have 
used ABFOW structure functions [@]. We see that for the low center-of-mass 
energy experiments WA70 gjj (VS = 23 GeV) and E706 @ (VS = 31 GeV) 
the gluon and the valence quarks contributions are of the same order at low 
Pu, whereas when P t i becomes larger, the valence quarks dominate. For ISR 
experiments |21|,|2^| (\fS = 63 GeV) the glue contribution dominates up to 



Pu — 10 GeV. For the UA2 experiment |23|], when the pseudo rapidity rj = 1.4, 
the glue contribution is important up to Pu — 35 GeV. Finally for LHC, in 
the Pu range between 30 and 1000 GeV the glue contribution represents (60 - 
80) % of the partonic cross-section. In all cases the "sea" contribution (s,c) is 
always negligible. 

In order to estimate the z range we are sensitive to we will study in table I 
the integrand of eq. (||), i.e.: 

< z >= J ; - V f> d dip < (is) 

fj?EiDf(z,M])Ei d -^ 
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with z varying between 2E 7T o/y/S and 1. Note that the partonic cross-sections 
reach their maximum for z = 1 while the fragmentation functions decrease with 
z. As we can infer from Table I the large z region is kinematically favored. We 
have used set I of fragmentation functions which will be discussed later. 

As already mentioned the fragmentation functions are known less accurately 
than the structure functions. Up to now they have been extracted at LO || from 
e + e~ annihilation and semi inclusive deep inelastic muon nucleon scattering. 
No extraction from hadronic colliders data has been performed so far. In the 
following we will carry out an extraction of tt° fragmentation functions at NLO 
accuracy, using three different approaches. 



3 Extraction of 7r° fragmentation functions. 

3.1 Selection of experimental data. 

We first discuss the experimental data we will use to extract the 7r° fragmen- 
tation functions. We first consider e + e~ collisions. The JADE collaboration 
\W has published data at y/S = 14, 22.5 and 34.4 GeV. We use the data at 



34.4 GeV, covering mainly the low zjj range (up to Zjj = 0.209). Data from 

the TPC collaboration 1171 at y/S = 29 GeV are given as — — Jr , therefore a 

' — v ° °~had paz H ' 

value of R=4.00 is assumed to bring them to the usual form ^^p- Data from 
the TASSO collaboration |l|] at yJ~S = 34.6 GeV extend up to z H = 0.728. The 
broadest zjj range is covered by data from the CELLO collaboration, extending 
from z H = 0.049 to z H = 0.919 at VS = 35 GeV g|] and from Z H = 0.094 to 
z H = 0.847 at VS = 22 GeV 0. Data from experiments at DORIS are not 
used, as hardly any point survive with the cut on the lower energy of the 7r° at 
2 GeV. Data obtained at LEP are for the moment not constraining. However, 
cross checks have been performed with the 2 points surviving the cut of data 
from the Argus collaboration [25| at \/~S = 10 GeV and the 4 points from the 
L3 collaboration p§ at \/S = 91 GeV. 



Let us discuss now experimental data from hadronic colliders. Data in 
hadronic reactions have been selected for this study taking into account sta- 
tistical and systematic accuracy. Whenever possible, reconstructed ir° are pre- 
ferred. For SPS fixed target energies, the available data in pp reactions are in 
reasonable agreement and we will use the data in the central rapidity range 
at \fS = 23 GeV, from the WA70 collaboration p|. The FNAL fixed target 
range overlaps with the lower ISR energy range. The recent data at \/S = 31 
GeV from pBe reactions obtained by the E706 collaboration [20] are in agree- 
ment with some of the ISR results. Resolved 7r° at y/S = 62.8 GeV taken from 
table 5 ( more precisely data corresponding to the super-retracted geometry) 
of Kourkoumelis et al. [g2| are used. They will be compared with other data 
available at this energy. We will use also the more recent data from the AFS 
collaboration [^], which however show a different Pt dependence. At collider 
energies, the latest data from the UA2 experiment at \/~S = 630 GeV with 
average pseudo rapidity rj = 1.4 will be used [23]. Cross checks have been made 
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with data at = 540 GeV with average pseudo rapidity i] = although ir° 
are not disentangled from direct photons. 

3.2 Fragmentation functions from HERWIG. 

We first consider the 7r° inclusive production in e + e~ annihilation at Mjq = 
\/~S = 30 GeV, as simulated by the Monte Carlo generator HERWIG. As well 
known, this event generator includes the QCD parton shower to leading and 
next to leading accuracy - in particular the kinematical corrections due to the 
phase space boundaries are summed up to all orders - as well as the hadro- 
nisation of the color singlet clusters into the physical particles. Furthermore 



HERWIG has been shown [24] to describe with good accuracy the observed fea- 
tures of PETRA and LEP data. Then we will use the 7r° distribution generated 
by each quark flavor which originates from the photonic vertex, as a realistic 
description of the quark fragmentation into ir°. Owing to the symmetry of 
quarks and antiquarks fragmenting into 7r° we extract the quark fragmentation 
functions from: 



(z H , M%) ~ 6a Y, 4 D l M /o)> ( 19 ) 



where the pointlike cross section <jq has been previously defined. The reaction 
e + e _ — > 7T + X has been therefore decomposed into each contribution e + e~ — ► 
uu, dd, ss, cc and bb. The generated distributions are parametrized as 

Df(z,M] ) = N lZ a >(l-zf> (20) 

and analyzed using the minimization procedure MINUIT. The coefficients iVj 
are constrained by the normalization condition: 

J^dzD l (z,M%) = (n 7T ) l , (21) 

where the average values {n- K ) i are given by HERWIG for each quark flavor, in 
agreement with the total observed multiplicity {n^). The parameters iVj, and 
Pi are extracted from the 7r° inclusive distribution generated, for each flavor, 
in the x range .025 < zh < -95 and shown in table II. As can be inferred from 
this table the statistical error on the parameters is less than 5%. 

As an illustration of the accuracy of the method and also of its limitations, 
the 7T° inclusive cross-section obtained from eqs. (0) and @, together with the 



results of table II , are compared in figure 2 with the CELLO data [15] at yS 
35 GeV. The agreement is reasonable in the range zh < -5. So far we have not 
included the contribution from the gluon fragmentation function. Indeed from 
the analysis of the three jet events it would be possible, in principle, to extract 
from HERWIG the appropriate information. The corresponding accuracy is 
however unsatisfactory, due to the limited sensitivity to hard gluon effects in 
e + e~ annihilation. 
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For this reason we have followed a different approach. To extract the gluon 
fragmentation function from HERWIG we have analyzed the subprocess gg — > 
gg — » 7T° + X from pp annihilation at Mjo = y/s ~ 30 GeV, in analogy to the 
quark case. In order to eliminate the background from the fragmentation of 
the spectator partons we have considered the pions lying only within a cone 
of semi aperture 5 = .35 — .40 rad around the direction of the parent gluons 
emitted at 90 deg. The value of 5 is found by an appropriate angular study of 
the generated distribution. With a parametrization of the form (^) we find the 
values of the parameters N g , a g and /3 g given in table III. After inclusion of the 
gluon fragmentation function and use of NLO evolved fragmentation functions 
together with NLO terms in the ir° inclusive cross section the agreement with 
CELLO data is improved as can be inferred from figure 3 up to zh — 0.8 



We compare now our predictions at NLO to experimental data from hadronic 
colliders. We first consider the data from CERN ISR ffi^ H , for y/S = 63 GeV, 



compared in figures 4 and 5 with our predictions for \i = M = Mt = Pt and 
[i = M = Mf = Pt/2 using the quark fragmentation functions from table II and 
the two gluon sets from table III, with 5 = 0.35 and 5 = 0.40. The agreement 
is satisfactory within the theoretical and experimental uncertainties. 



Let us focus now on the UA2 data at the SppS collider |23|]. We will use 
two sets of quite precise data, for P t < 15 GeV and r\ ~ and, for 15 < Pt < 45 
GeV and rj ~ 1.4. The comparison with the theoretical predictions is shown 
in figures 6 and 7 for fi = M = Mf = Pt/2, Pt and for the two gluon sets 
of fragmentation functions. The agreement is quite good, and slightly favors 
the set corresponding to 5 = 0.35. The dependence on the renormalization, 
factorization and fragmentation scales at NLO will be discussed later. 



In the next subsection we will extract the ir° fragmentation functions at 

/o 



next to leading order using two different hypotheses at the reference scale M 2 



2 GeV 



3.3 Set I: fragmentation functions with natural scales. 

For this set, we take a s as given by equation (|i~4|) and A = 190 MeV, 
corresponding to the set of structure functions we will use |[ . 



3.3.1 Definition 

We assume for this case an SU(2) symmetry: 

DZ\z,M%) = D£(z,M%) = Df{z,M%) = D£(z,M%) = D v {z,M%)+D s {z,M%) 

(22) 

Then, we take 



Df(z,M%) = Df(z,M%) = Df(z,M%) = D?(z,M%) = D s (z,M%), 

' (23) 
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and 



Dg (z,Mf ) = D G (z,Mf ). (24) 
We parametrize the different functions of z as follows 

D v (z,M%) = N v (l-zf« (25) 
D s {z,M%) = N s (l-z)P° (26) 
D G {z,M%) = N g (1 — z)^ 9 . (27) 

At the initial scale Mjq, we start with four flavors. The b quark contribution is 
taken into account in the evolution. Fixing the threshold at 4 m 2 , so we have: 

D b (z, M f ) - | ^ (1 _ z)Pa . f M i = 4 ^ (28) 

So we are left with six parameters to be determined with the help of experi- 
mental data. 



3.3.2 Choice of the scale 

We use the standard approach to fix all the scales to the same value which 
is some natural scale of the problem. More precisely, for e + e~ collisions, we 
take \i = Mf = y/S whereas for p p collisions, we set the three scales equal and 
proportional to the transverse momentum of the 7r° : 

H = M = M } = cP t 

where c is a constant to be fixed by the fit to experimental data . 



3.3.3 Results for set I 

First of all, for e + e~ collisions, we limit ourselves to a ir° energy greater 
than 2 GeV because we don't trust perturbation theory for low 7r° energies. 
Therefore for y/S ~ 30 GeV, we will only use z values greater than 0.1. As 
it can be inferred from eqs (25-27) we have not used a factor z a in the input 
parametrizations since in this z range it does not improve the fit but only leads 
to correlations. With six parameters, a big correlation still occurs between N v 
and (3 V , so we fix (3 V = 1. Then N s , N g and j3 g remain slightly correlated. A 



good fit to CELLO @, TASSO ||], TPC and JADE [gj] data leading to 
a x 2 = 26.3 for 29 points is obtained for values of the parameters given in table 
IV (systematic errors have been added in quadrature to statistical errors). 



Using set I of fragmentation functions we will now evaluate the NLO cross- 
sections for inclusive ir° production in hadronic collisions and compare them to 
experimental data from low center of mass energies up to the CERN collider 
one. Here, the situation is less clear. First, if we keep constant the value 
of the parameter c it is impossible to obtain a good fit in the whole energy 
domain. For example, setting c ~ 1.5, the ISR data can be described but the 
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theoretical predictions are by far too low for WA70 and E706 and too high 
for UA2. A simple solution to this problem is to allow c to vary with the 
hadronic kinematical variables, in particular y/S. A correct description of the 
data requires c ~ 0.39 for WA70 |J] (see figure 8a), c ~ 0.5 for E706 ||](see 
figure 8b), c ~ 1.5 for ISR experiment pT| , p2| (see figure 8c) and c ~ 5.5 for 
UA2 |^3[ (see figure 8d). In particular for the ISR energy range, the data from 
AFS collaboration |2lj] are marginally consistent with those of reference p2~| 
since the transverse momentum dependence in the two experiments is different. 
Therefore it is very difficult to describe both ISR data with high precision. We 



get rather good fits of data of Kourkoumelis et al. [22] with x = 20.6 for 
14 points using /j, = M = Mf = l.SPt and of the AFS collaboration |2l]] with 



X 2 = 12.2 for 11 points using \x = M = Mr = 1.6Pf. Notice that the slope of the 



UA2 data is not correctly reproduced, with a x 2 = 50.2 for 11 points. The x 2 
have been calculated with statistical errors, allowing the overall normalization 
to vary within the systematic error. 

A comment is in order here. The approach followed so far is rather simple. 
When the energy grows up the scales needed to describe data have also to 
increase. As stated above an acceptable fit of UA2 data [23] in the forward 
direction can be obtained for the choice of scales [i = M = Mf = 5.5Pt which 
is a priori a large scale. The compensation occurring between the leading and 
next-to-leading terms concerning the scale dependence is much more effective 
at high energies. At low energy, since we prevent the scale to be less than 
Mf Q = y/2 GeV, this compensation does not occur and the behavior of the 
leading and next-to-leading cross-sections is quite the same. In other words, we 
are not in a good region to perform perturbation theory. 

This approach might be criticized. Indeed it is not very predictive, since 
the scales change with the energy. In other words one adds a new parameter 
which acts as an overall normalization for each experiment. Notice that the 
normalization of the glue fragmentation function N g is strongly correlated to 
the choice made for the scale. More precisely, we could perfectly find a value 
for N g which describes the UA2 data with c = 0.5. But in this case we couldn't 
describe the other data at lower energies. 



3.4 Set II: fragmentation functions with optimized scales. 

For this set, we take the numerical solution of equation (|1|) for a s and 
A = 230 MeV, since we will use the ABFOW set of structure functions . 

3.4.1 Definition 

We assume also for this case an SU(2) symmetry: 

Dl\z,M%) = Df(z,M 2 f0 ) = Df(z,M 2 fQ ) = Df(z,M 2 ) = D u (z,M 2 ). 

(29) 

Then we take: 

Df(z,M 2 f0 ) = Df(z,M 2 ) = D s (z,M 2 f0 ), (30) 
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Dl (z,M 2 f0 )=D? (z,M 2 f0 ) = D c (z,M 2 f0 ), (31) 

and 

Df(z,M%) = D g (z,M%). (32) 
We parametrize these different functions of z in the following way: 

(33) 
(34) 
(35) 
(36) 

So we are left with eight parameters to be determined with the help of experi- 
mental data. Since we will use the optimized procedure for the determination 
of the scales, it is much simpler not to change the number of flavors. So, in this 
case, we will neglect the b contribution. This assumption is motivated by the 
fact that a(e + e~ —* 7* —* b b) = 1/4 cr(e + e~ — > 7* —* c c) and in pp collision 
the b production is suppressed due to the weak b content of the proton. 

A few remarks are in order here. As in the case of set I, the non singlet 
part D~ is always zero due to our assumptions. We did not take = D£ 
because in this case the sum over the four flavors of Df weighted by the square 
electric charge is zero: 



D u (z, 


M%) = N u 


z~ 


1 (1 


-zf u 


D s (z, 


M%) = N s 


z~ 


1 (1 


-zf> 


D c (z, 


M%) = N c 


z~ 




-z)^ 




M%) = N g 


z~ 




-zf» 



E 



e 2 (Df(z,M 2 ) + D+(z,M 2 ))=0. 



So, there is no non-singlet contribution to the cross-section. Therefore we could 
parametrize directly the singlet and the glue with four parameters only. The 
e + e~ data could be correctly described, but the glue is very constrained and it 
will not be possible to fit hadronic data in the whole energy range. 

3.4.2 Choice of the scale 

For set II, we use optimized scales according to the procedure of Politzer 
and Stevenson [28|. Concerning e + e~ collisions, our approach is the following. 



Firstly since the scale does not appear at lowest order, we cannot optimize 
with respect to it. Therefore we set \i = Mf and perform an optimization only 
with respect to the scale Mf. Therefore, a priori, our optimized scale depends 
on the choice made for the input fragmentation functions. We have not found 
a way to get rid from this sensitivity. In practice, the optimized point changes 
slowly when the input is modified and in addition, since we are in a stable 
region, it does not matter if we are not exactly on the optimized point. The 
optimized scale Mj pt is of order of \/~S/8 varying slowly with z. Furthermore, 
we find no optimization scale for z < .03 for y/S = 35 GeV, z < .05 for y/~S = 29 
GeV and z < .1 for yS = 22 GeV. For lower values of yS, it is not possible to 
optimize. 
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We also use an optimization procedure for hadronic collisions. So we require 
that: 

= (37) 

d E /_^±2^1 = o (38) 



91n(MVA2) * tfp^ 

E /_5l±2^1 = o. (39 ) 



91n(M|/A2) - ^ 
The first equation can be computed analytically: 

e ^p+p^ = _aV)& I ^ A + 3(1 + b'a s (n 2 )) 



ain(/xVA 2 ) w ^ 



2 M ln(g 



(40) 



having used 



mn(l , /x2) =-ba 2 s (» 2 )(l + b'a s (» 2 )). (41) 

Note that terms of order of have been cancelled as it should be. Now, we 
determine the scale \i in order to cancel the right-hand side of eq (|40|). This 
ensures us that the corrective term K will be negative with a magnitude of 
roughly 10 % of the lowest order. Then we compute numerically the value of 
the scales M and Mj which have to fulfill the equations ( |3l| ) (|39|), the scale \i 
being now a function of M, Mf. We require that the factorization scales must 
be greater than v2 GeV and that the renormalization scale is such that the 
running coupling constant a s is less than .34. With these constraints it will 
be impossible to optimize in low Pt range. More precisely, for low center of 
mass energies (y/S < 63 GeV), the optimization is not possible for Pt < 5 GeV. 
Therefore these regions are not appropriate to apply an optimization procedure. 

3.4.3 Results for set II 

First we freeze /3 S , (3 C and (5 g according to the counting rules. There are still 
too many parameters, so we fix N g and fit to e + e~ data with four parameters 
N u , (3 U , N s and N c . The fragmentation functions extracted are then used to 
evaluate hadronic cross sections. Then we vary N g refitting e + e~ data and 
apply the new input to pp data. This procedure is repeated until a reasonable 



description of hadronic data is reached. Good fits of e + e data (CELLO [1J, 15|, 



TASSO P|, TPC and JADE @) leading to a X 2 ^ 1 per d.o.f. are 



obtained for the two sets - hereafter denoted as set Ha and set lib - displayed 
in Table V and Table VI (see figures 9 and 10 using set lib). The two sets differ 
mainly for the gluon normalization. As can be seen from inspection of figures 
11a, lib, 12a and 12b a rather good fit of the latest UA2 data at y/S = 630 



GeV |23J, AFS J21J] and Kourkoumelis et al data []22j can be obtained leading 



to a x — 50 for 31 points. Kourkoumelis et al. data favor the set characterized 
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by the largest glue (set lib) whereas UA2 data are better fitted by the other 
set (set Ha). Notice that we have taken into account the systematic errors of 
the data which affect the overall normalization. The x 2 are 3.46 (4.28) for the 
11 AFS points, 31.54 (23.52) for the 9 Kourkoumelis et al. points and 14.91 
(20.00) for the 11 UA2 points with the parameters of set Ha (lib). Inside the 
systematic errors we can also describe UA2 data at yS = 540 GeV and rj = 1.4. 
On the other hand we are not able to describe WA70 and E706 data with the 
values of N g found before. This is not very surprising since the corrective term 
is found to be huge, and although we can find an optimization point this is 
not very stable suggesting that we are not in the appropriate region to trust 
perturbation theory. 

4 Predictions at future hadron colliders. 

As we have seen present data do not allow to extract the ir° fragmenta- 
tion functions unequivocally. To this aim the forthcoming information from 
ep HERA collider should be very helpful. With these limitations we will now 
estimate the 7r° rates at LHC using the various sets of fragmentation functions 
previously derived. 

Let us consider first set I of fragmentation functions. In order to describe 
hadronic data we had to increase the scales \i = M = Mf from -y at \/~S ~ 20 
GeV up to 5Pt at ^/S = 630 GeV. An extrapolation to LHC energy would lead 
to fx = M = Mf ~ 50Pt which seems by far an unnatural scale. To estimate 
the sensitivity to scales we show in figure 13 the ratio of cross sections at LHC 
for the two scales 50Pt and Pt at rj = 0. As can be inferred from the figure 
the rates differ by at most a factor of two. To estimate the uncertainty due to 
structure functions we have taken the set of structure functions of HMRS || 
using the MS scheme and the set of Morfing-Tung Q using the DIS scheme. 
The predictions differ by at most 20%. Similarly the ratio of predictions using 
set II is displayed in figure 14. 

Finally the ratio of the two predictions from the HERWIG fragmentation 
functions, for 5 =0.35, 0.40, evolved to NLO accuracy as discussed in section 
3.2, are displayed in figure 15. 

The situation is summarized in figure 16 where we show the absolute rates 
at LHC for ri = from the most plausible sets in the three approaches. This 
gives an estimate of the theoretical uncertainty which is of the order of a factor 
two. The uncertainty on structure functions is marginal compared to the poor 
determination of fragmentation functions. 

To show the stability of the NLO corrections we display the cross section as 
a function of the scales \i and M = Mf compared to the LO result for Pt = 50 
GeV (figures 17) and for P t = 200 GeV (figures 18). We vary the scales between 
Pt/5 and 5Pt. The NLO cross sections exhibit a saddle point whereas the LO 
cross sections decrease monotonically when the scales increase. 
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The uncertainty due to factorisation scheme, especially coming from frag- 
mentation functions is expected to be tiny for the two following reasons. Firstly 
the evaluation done for one jet inclusive cross section has shown [2] that at col- 
lider energies its magnitude is of the order of 5% -if done correctly - and we 
can reasonably expect a same order of magnitude for one hadron inclusive cross 
section. Secondly a precise estimate doesn't seem mandatory compared to the 
large uncertainty coming from fragmentation functions. 
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5 Conclusions 



We have performed a complete next to leading order analysis of inclusive 
7T° production from e + e~ and hadronic data. For the first time an attempt to 
extract sets of 7r° fragmentation functions at NLO has been performed. The 
present quality of data does not allow us to derive a unique set fitting all ex- 
perimental data. For this purpose more accurate measurements from hadronic 
colliders in the large Pt domain, from e + e _ colliders in the large Zjj domain as 
well as complementary information from ep collisions will be very helpful. The 
theoretical uncertainties are mainly due to the poorly determined fragmenta- 
tion functions. Nevertheless the absolute rates at future colliders like LHC and 
SSC can be predicted within a factor of two. This will certainly be of help for 
neutral background rejection at supercolliders. 



Note added in proof 

After completion of this work the paper "Higher order QCD corrections to 
inclusive particle production in pp collision" by F. M. Borzumati et al. [^] has 
appeared, where the 7r° inclusive production has been discussed, using the old 
LO fragmentation functions of ref.[8] and the NLO results of our group 0. 
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Table captions 



• Table I: average fraction of energy of the fragmenting parton (see eq.18) 
from pp collisions. 

• Table II: parameters of the quark fragmentation functions (see eq.20) as 
obtained from HERWIG in e + e~ annihilation at M = 30 GeV. 

• Table III: parameters of the gluon fragmentation functions (see eq.20) as 
obtained from HERWIG at Mq = 30 GeV, with two hypotheses on the 
angle S (see text). 

• Table IV: parameters for the 7r° fragmentation functions (set I) obtained 
from e + e~ annihilation. 

• Table V: parameters for the 7r° fragmentation functions (set Ha) obtained 
from optimization procedure in e + e _ annihilation. 

• Table VI. Parameters for the 7r° fragmentation functions (set lib) obtained 
from optimization procedure in e + e~ annihilation. 
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y/S = 23 GeV and n = 0. 


y/S = 63 GeV and n = 0. 


y/S = 630 GeV and n = 1.4 




< z > 




< z > 




< z > 


4.11 


0.81 


5.25 


0.67 


13 


0.55 


4.61 


0.82 


6.73 


0.70 


21 


0.60 


5.69 


0.86 


8.23 


0.73 


29.8 


0.65 


6.69 


0.89 


10.4 


0.77 


43.7 


0.74 



Table I 



Process 


a 


(3 




< n w > 


e + e~ — ► uii 


-0.95 ±0.02 


3.67 ±0.19 


1.20 


2.95 


e + e~ — > cid 


-0.95 ±0.02 


3.67 ±0.15 


1.24 


2.87 


e + e~ — > ss 


-0.88 ± 0.02 


5.32 ±0.23 


1.68 


2.73 


e + e~ — > cc 


-0.82 ± 0.02 


8.02 ±0.24 


3.09 


3.42 


e + e~ — > 66 


-0.95 ±0.02 


10.94 ±0.29 


2.92 


4.20 



Table II 



5 


a 


(3 


N g 


< n w > 


0.35 rad 


-0.28 ±0.04 


6.71 ±0.39 


14.49 


3.65 


0.4 rad 


-0.37 ±0.04 


5.79 ± 0.36 


12.93 


4.55 



Table III 



Parton 


Cii 


ft 




valence 


0. 


1. 


0.19 


sea 


0. 


5.2 


3.5 


gluon 


0. 


2.03 


4.9 



Table IV 



Parton 


Cti 


Pi 




up 


-1. 


0.94 


0.11 


strange 


-1. 


3.0 


0.55 


charm 


-1. 


4. 


2.7 


gluon 


-1. 


2. 


0.55 



Table V 
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Parton 




ft 




up 


-1. 


1.11 


0.15 


strange 


-1. 


3.0 


0.18 


charm 


-1. 


4. 


2.5 


gluon 


-1. 


2. 


0.75 



Table VI 
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Figure Captions 



• Fig. 1: relative inclusive partonic production in hadronic collisions as 
a function of the partonic transverse momentum P t \ at various ener- 
gies: WA70 experiment (fig la), ISR experiments(fig lb),UA2(fig lc) and 
LHC(fig Id). The various curves refer to: pp — > u+d (full line), pp — > s + c 
(dashed line) and pp — ► 5 (dot-dashed line) . 

• Fig. 2: LO inclusive 7T° production in e + e~ annihilation with the quark 
fragmentation functions extracted from HERWIG at \/S = 30 GeV, com- 
pared with experimental data at V~S = 35 GeV. 

• Fig. 3: NLO inclusive 7r° production in e + e~ annihilation with the quark 
and gluon fragmentation functions extracted from HERWIG and evolved 
at \/S = 35 GeV, compared with data. 

• Figs. 4: NLO inclusive 7r° production in pp collisions at ISR energies for 
jj, = M = Mf = P t ,P t /2 for 5 = 0.35 (see text). The fragmentation 
functions are from HERWIG. 

• Figs. 5: same as fig 4, for 5 = 0.40. 

• Figs. 6: NLO inclusive 7r° production in pp collisions at SppS energies 
for n = M = M f = P t , P t /2 for 6 = 0.35 and S = 0.40,at y/S = 540 GeV 
and 77 = 0. The fragmentation functions are from HERWIG. 

• Figs. 7: same as fig 6 at y/S = 630 GeV and 77 = 1.4. 

• Figs. 8: NLO inclusive 7r° production in hadronic collisions with set I 
fragmentation functions (see text) for various energies. The scales /_t = 
M = Mf = cPt are indicated explicitely. 

• Fig. 9: NLO inclusive 7r° production in e + e~ annihilation with set II of 
fragmentation functions compared to CELLO, TASSO and JADE data 
. The gluon parameter N g takes the value N g = 0.75. Data and theory 
have been multiplied by 0.1 at x/~S = 22 GeV. 

• Fig. 10: Same as fig 9 for TPC data. 

• Figs. 11: NLO inclusive tt° production in pp collisions with set II frag- 
mentation functions (see text) at ISR energies (squares correspond to AFS 
data whereas circles correspond to Kourkoumelis et al data). The gluon 
parameter N g takes the value N g = 0.55 (fig 11a) and N g = 0.75 (fig lib). 

• Figs. 12: same as fig 11 at SppS energies. 

• Fig. 13: ratio of inclusive 7r° cross sections predicted at LHC for i] = 
using set I of fragmentation functions with fi = M = Mf = 50Pt over 
H = M = M f = P t . 



19 



Fig. 14: ratio of inclusive n cross sections predicted at LHC for r] = 
using set lib of fragmentation functions over set Ha . 

Figs. 15: ratio of inclusive 7r° cross sections predicted at LHC for i] = 
using fragmentation functions from HERWIG. 

Fig. 16: inclusive 7r° cross sections predicted at LHC for r/ = using 
various fragmentation functions : HERWIG with 5 = 0.35 (full line), set 
I with fj, = M = M f = 50P t (dot-dashed curve) and set II with N g = 0.75 
( dashed curve). 

Figs. 17: inclusive 7r° cross sections at LHC in pb using set I of fragmen- 
tation functions as a function of the scales \i and M = Mf for P t = 50 
GeV and rj = 0. LO prediction : fig 17a . NLO prediction : fig 17b. 

Figs. 18: same as fig 17 for P t = 200 GeV and rj = 0. 
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